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We give a theoretical analysis of published experimental studies of the effects of impuri- 
ties and disorder on the superconducting transition temperature, T c , of the organic molec- 
ular crystals k- (BEDT-TTF) 2 X (where X=Cu[N(CN) 2 ]Br and Cu(NCS) 2 and BEDT-TTF is 
bis(ethylenedithio)tetrathiafulvalene) and /3-(BEDT-TTF) 2 X (for X=I 3 and IBr 2 ). The Abrikosov- 
Gorkov (AG) formula describes the suppression of T c both by magnetic impurities in singlet su- 
perconductors, including s-wave superconductors and by non-magnetic impurities in a non-s-wave 
superconductor. We show that various sources of disorder (alloying anions, fast electron irradiation, 
disorder accidentally produced during fabrication and cooling rate induced disorder) lead to the 
suppression of T c as described by the AG formula. This is confirmed by the excellent fit to the 
data, the fact that these materials are in the clean limit and the excellent agreement between the 
value of the interlayer hopping integral, t±, calculated from this fit and the value of t± found from 
angular-dependant magnetoresistance and quantum oscillation experiments. There are only two 
scenarios consistent with the current state of experimental knowledge. If the disorder induced by all 
of the four methods considered in this paper is, as seems most likely, non-magnetic then the pairing 
state cannot be s-wave. We show that published measurements of the cooling rate dependence of 
the magnetisation are inconsistent with paramagnetic impurities. Triplet pairing is ruled out by 
NMR and upper critical field experiments. Thus if the disorder is non-magnetic then this implies 
that I > 2, in which case Occam's razor suggests that d-wave pairing is realised in both /3-(BEDT- 
TTF) 2 X and k-(BEDT-TTF) 2 X. However, particularly given the proximity of these materials to 
an antiferromagnetic Mott transition, it is possible that the disorder leads to the formation of local 
magnetic moments via some novel mechanism. Thus we conclude that either /3- (BEDT-TTF) 2 X 
and k-(BEDT-TTF) 2 X are d-wave superconductors or else they display a novel mechanism for the 
formation of localised moments, possibly related to the competition between the antiferromagnetic 
and superconducting grounds states. We suggest systematic experiments to differentiate between 
these two scenarios. 

PACS numbers: 74.20.Rp, 74.62.-c, 74.70.Kn 



I. INTRODUCTION 

Superconductivity is often found near magnetic order- 
ing. This may be antiferromagnetic (AFM) order such as 
in the cupratesi and the heavy fermion superconductors^ 
or ferromagnetic order as in the ZrZii2 or UGe2 (see Refs. 
0and0 respectively). In each of these cases it is believed 
that the superconductivity is unconventional ^LLSi that 
is to say that the Cooper pairs have a non-zero angu- 
lar momentum. The issue of unconventional supercon- 
ductivity near magnetic ordering is of general interest 
because it may lead to insights into both non-phononic 
pairing mechanisms- and the theory of quantum critical 
points^ 

Despite the fact that it is now twenty years since su- 
perconductivity was discoverediiii^ in the layered or- 
ganic compounds (BEDT-TTF) 2 X (where BEDT-TTF is 
bis(ethylenedithio)tetrathiafulvalenc and X is an anion, 
e.g., Cu[N(CN)2]Br or I3) the pairing symmetry remains 
a matter of debated BEDT-TTF salts form a num- 
ber of crystal structures which are denoted by a Greek 
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letter. All of the crystal structures consist of alternat- 
ing layers of BEDT-TTF and an anionic In /3-(BEDT- 
TTF) 2 X and k-(BEDT-TTF) 2 X, which we consider here, 
the BEDT-TTF molecules form a dimerised structure 
where the anion removes one electron per dimer. Thus 
we have alternating conducting (BEDT-TTF) and insu- 
lating (anion) layers. A particularly interesting feature 
of these materials is that they can be driven from an 
AFM insulating state to a superconducting state by the 
application of hydrostatic pressure or by changing the 
anion. 15,16 

In principle the simplest way to identify the pair- 
ing symmetry, or at least the nodal structure, of a su- 
perconductor is to measure the low temperature be- 
haviour of thermodynamic or transport properties. For 
example the specific heat follows an exponentially acti- 
vated temperature dependence for a nodeless gap (Cy oc 
exp(— |A(0)|/fcsT), where A(0) is the superconducting 
gap at zero temperature) and a power law dependence for 
a gap with nodes (CV cx T 2 for line nodes and CV cx T 3 
for point nodes on a three dimensional Fermi surface). 17 
In practice however, there are difficulties associated with 
this method of identifying the pairing symmetry, not the 
least of which is the need to make measurements at ex- 
tremely low temperatures. (Typically a wide tempera- 



ture range is required in the region T/T c < 0.2, so in 
the case of /c-(BEDT-TTF) 2 Cu[N(CN) 2 ]Br (T c - 10 K) 
one requires measurements taken over a wide range of 
temperatures below ~ 2 K.) The apparently strong 
couplingiSiiS nature of the superconductivity in theses 
charge transfer salts means that the behaviour of ther- 
modynamic and transport functions near T c is unable to 
differentiate between pairing states on symmetry grounds 
alone and so we must wait for calculations based on a 
specific theory of superconductivity to use this data to 
examine the pairing symmetry. 

Regardless of the reasons one fact is clear ^ low tem- 
perature behaviours have been, to date, unable to settle 
the debate on the pairing symmetry in the layered or- 
ganic superconductors. In particular, two pairing sym- 
metries have been widely discussed: strong coupling s- 
wave superconductivity and d-wave pairing. 

In K-(BEDT-TTF) 2 Cu[N(CN) 2 ]Br the 13 C NMR spin 
lattice relaxation rate)22i2ii2£ (Ti) _1 , shows no Hebel- 
Slichter peak and a power law cutoff, (Ti) _1 oc T n , where 
n ~ 3. A Hebel-Slichter peak is expected for s-wave 
pairing while (Ti) _1 oc T 3 is expected for line nodesi^ 3 . 

Much controversy has surrounded the London pen- 
etration depth with some groups reporting s-wave 
pairingS^S^LSLS^ and others finding line nodes con- 
sistent with d-wave pairing 29 ; 30 ! 31 ! 32 ! 33 ^ in both k- 
(BEDT-TTF) 2 Cu[N(CN) 2 ]Br (Refs. I24l25l2fil and || 
pl3pp) and ^(BEDT-TTF) 2 Cu(NCS) 2 (Refs. I2l2rl 
1271281 r i3lL I30L I33I and 13^) . However, the most recent 
measurements^! have two advantages over older experi- 
ments. Firstly very low magnetic fields were used. The 
use of fields less than the lower critical field is important 
in penetration measurements because vortex dynamics 
are a serious impediment to accurately measuring the 
penetration depth. Secondly Carrington et alm^ made 
measurements down to 0.4 K and therefore made a large 
range of measurements below T ~ 0.2T C . This is the 
lowest temperature range considered in any of the ther- 
modynamic or transport experiments, making Carring- 
ton et aZ.'s conclusions the most reliable drawn from ex- 
periments of this type. Carrington et al. found that 
the temperature dependence of the penetration depth 
of both K-(BEDT-TTF) 2 Cu[N(CN) 2 ]Br and k-(BEDT- 
TTF) 2 Cu(NCS) 2 is inconsistent with a nodeless gap. 

Initial measurements of the specific heat of k-(BEDT- 
TTF) 2 Cu[N(CN) 2 ]Br showed a T 2 dependence 3 ^ but the 
interpretation of these results has been questioned^ 
More recent measurements of the specific heat have found 
an exponentially activated temperature dependence for 
both K-(BEDT-TTF) 2 Cu[N(CN]2jBr (Ref. HI) and n- 
(BEDT-TTF) 2 Cu(NCS) 2 (Ref. M). 

Several groups have considered probes which do not 
rely on the low temperature behaviour of the measure- 
ment. Brando et til.— and Arai et a/i^LS attempted to 
observe the local density of states (LDOS) of k-(BEDT- 
TTF) 2 Cu(NCS) 2 by measuring the differential conduc- 
tance using a scanning tunnelling microscope (STM). 
Each of these experiments found a LDOS that is con- 



sistent with d-wave pairing however, none of the ex- 
periments observed the coherence peaks which are a 
characteristic feature of the superconducting state and 
have been observed^ifiiii in similar experiments on 
Bi 2 Sr 2 CaCu 2 08+a;. Also one should note that Bando 
et al— observed a LDOS in the layered s-wave super- 
conductor NbN which has the same form as that which 
is interpreted as d-wave in experiments on k-(BEDT- 
TTF) 2 Cu(NCS) 2 . 

Schrama et al— attempted to determine the 
anisotropy in the superconducting order parameter by 
measuring the magneto-optical properties of k-(BEDT- 
TTF) 2 Cu(NCS) 2 and found results indicative of d-wave 
pairing. However, in light of the debate over the inter- 
pretation of these results^ 3 ^^ 5 . one cannot consider these 
measurements to have determined the pairing symmetry. 

Izawa et alm^ measured the thermal conductivity ten- 
sor of k-(BEDT-TTF) 2 Cu(NCS) 2 in a magnetic field. 
They observed a four fold anisotropy at low tempera- 
tures which they interpreted as evidence for d-wave pair- 
ing. However, it is possible that the vortices produced in 
k-(BEDT-TTF) 2 Cu(NCS) 2 are actually Josephson vor- 
tices. Therefore it remains to be shown whether or not 
the theory^L^Sii 9 . on which Izawa et al. base their analy- 
sis is valid for this material. 

The 13 C NMR Knight shift has been measured 2 ^ 
for K-(BEDT-TTF) 2 Cu[N(CN) 2 ]Br. With a magnetic 
field, H, parallel to the conducting planes, as T — > 
so does the Knight shift i* 3 - This does not actually rule 
out triplet pairing, although it does make triplet pair- 
ing extremely unlikely. This experiment is compatible 
with a triplet state in which d(k) xH = where d(k) is 
the usual Balian-Werthamer order parameter for triplet 
superconductivityi^Siiii An example of a triplet phase 
compatibles^ with this experiment is an A-phase with 
d(k) pinned to the c-axis^ 7 ^ which is not an impossi- 
bility given the highly anisotropic nature of k-(BEDT- 
TTF) 2 Cu[N(CN) 2 ]Br. However, Zuo et alM measured 
the critical field as a function of temperature with H 
parallel to the conducting planes. In this configuration 
no orbital currents flow so the critical field is due to 
Clogston-Chandrasekhar (or Pauli) limiti 53 i 55 i 56 There 
is no Clogston-Chandrasekhar limit for H _L c for triplet 
states compatible with measured Knight shift. Thus for 
such states there would be no critical field with H||6 
(in fact for such states one would increase T c by apply- 
ing a field parallel to the b-axis^ 3 -). Experimentally 57 
it is found that superconductivity is destroyed by a 
magnetic field parallel to the b-axis. Therefore only 
when considered together do the three experiments dis- 
cussed aboveSSi 5 ^ 7 . strictly rule out triplet pairing. 175 
Further evidence for Clogston-Chandrasekhar limiting 
comes from the observation that the in plane upper crit- 
ical field is independent of the field direction^ Given 
the anisotropic nature of the Fermi surface of k-(BEDT- 
TTF) 2 Cu[N(CN) 2 ]Br it is extremely unlikely that or- 
bital mechanisms for the destruction of superconductiv- 
ity would be so isotropic. 
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The results of quantum chemistry calculations sug- 
gest that the simplest theoretical model which can de- 
scribe these materials is a half-filled Hubbard model 
on an anisotropic triangular lattice>iSiS Because of the 
proximity of the antiferromagnetic insulating phase and 
the superconducting phase several groups have exam- 
ined the possibility of spin fluctuation induced super- 
conductivity within the confines of this model using 
a variety of techniques, including mean-field theory^ 
the fluctuation-exchange approximationj£fli£ii£2 third or- 
der perturbation theory^ weak coupling renormalisation 
group analysis pi the random phase approximation^^ 
and quantum Monte Carlo methods** 7 . All of these groups 
concluded that spin fluctuations lead to d-wave pairing. 
These authors found an enhanced dynamical suscepti- 
bility at (tt, ±7r) which leads to d x 2_y2 pairing. Alter- 
natively, both d-wave££ and s-waveSStBili pairing sym- 
metries have been considered in the context of phononic 
pairing mechanisms. 

So, perhaps the only emerging consensus is that the 
low temperature behaviours have not been able to con- 
clusively settle the debate between s-wave and d-wave 
pairing symmetries. In the remainder of this paper we 
will investigate how the effects of disorder can be used to 
distinguish between these two symmetries. 



II. THE ABRIKOSOV-GORKOV FORMULA 

Anderson's theorem 7 ^ states that for s-wave pairing 
non-magnetic impurities do not change T c . This is be- 
cause Cooper pairs are formed from time reversed states 
and although non-magnetic impurities may change, for 
example, the phonon spectrum, they do not break time 
reversal symmetry (TRS). However, magnetic impurities 
strongly reduce T c for all singlet states because they 
do break TRS. 73 This behaviour is described by the 
Abrikosov-Gorkov (hereafter AG) formula^ 



1 T c M2 4nk B T c T M 



(1) 



where T c q is the superconducting critical temperature in 
the pure system and ip{x} is the digamma function, tm is 
the quasiparticlc lifetime due to scattering from magnetic 
impurities. Assuming isotropic scattering tm is given 
byS 



tm 



= N M TrJi(Ji + l)N(0)\u M \ 



(2) 



where Nm is the number density of magnetic impurities, 
A(0) is the density of states per spin at the Fermi level, 
Ji is the total angular momentum of the paramagnetic 
atoms and um is the amplitude for scattering from a 
magnetic impurity. 

In the superconducting state the anomalous Green's 
function, F a p(k^u! n ), is finite and therefore there is, in 
the presence of non-magnetic impurities, an anomalous 



self energy, £2 afli^n), which, in n dimensions, is given 
byii 



1 



d n k 



2ttN(0)t n J (2tt) 



F a/3 (k,w„) (3) 



where r^r, the lifetime for scattering from non-magnetic 
impurities is given by2^ 



T/V 
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where Nn is the number density of non-magnetic im- 
purities and un is the amplitude for scattering from a 
non-magnetic impurity. 

For s-wave pairing £2,0/3(^,1) is clearly finite, and it can 
be shown that the anomalous self energy cancels exactly 
with the normal self-energy, £i ia/ g(w n ), when the criti- 
cal temperature is evaluated. Therefore T c is unchanged 
by non-magnetic impurities for an s-wave superconduc- 
tor, as expected from Anderson's theorem^ However, 
for non-s-wave pairingi£& it can be seen, from symmetry 
grounds alone, that the integral in equation © vanishes. 
Thus the anomalous self-energy does not cancel the nor- 
mal self energy and T c is lowered by non-magnetic impu- 
rities in a non-s-wave superconductor. Further, it can be 
shown that for pairing states with non-s-wave symmetry 
non-magnetic impurities reduce T c via the Abrikosov- 
Gorkov formulaii2i2i However, in this case 



(•5) 



where again we have assumed isotropic scatter- 
ing. The predictions 177 of Anderson's theorem 
have been confirmed for the alloys of many s-wave 
superconductors .iLZgjSLSI 

Hasegawa and Fukuyama-i suggested that weak local- 
isation could lead to an alternative mechanism for the 
suppression of T c in organic superconductors. Notably 
this mechanism allows for the suppression of T c by non- 
magnetic disorder in s-wave superconductors, in viola- 
tion of Anderson's theorem. However the Hasegawa- 
Fukuyama mechanism has a dramatically different r de- 
pendence to the AG formula. We will show in this paper 
that the observed suppression of T c in /3-(BEDT-TTF) 2 X 
and k-(BEDT-TTF) 2 X is described by the AG formula 
and therefore the predictions of Hasegawa and Fukuyama 
are not in agreement with experiment. For a multiband 
superconductor interband scattering processes can also 
lead to a suppression in T c (see, for example Ref. I82T1 . 
However, of the two polymorphs discussed in this paper 
only one (k-(BEDT-TTF) 2 X) has multiple sheets to its 
Fermi surface. As it seems reasonable to assume (unless 
evidence is found to the contrary) that the suppression 
of T c in both materials is due to the same mechanism 
we will not discuss interband scattering effects further. 
Also not that for moderate amounts of disorder interband 
scattering effects and the AG formula give very different 
predictions for the suppression of T c . 
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It can be shown that the digamma function has the 
property 



(6) 



Hence for H/t <C ksT c (i.e. as the number of impurities 
tends to zero) the AG equation becomes 



■kK l 



(7) 



Clearly the above is valid for both magnetic impurities 
in singlet states (r = tm) and non-magnetic impurities 
in non-s-wave pairing states (in which case r = r^r). 



A. Mixed order parameters 

As well as s-wave pairing and non-s-wave pairing, a 
third logical possibility exists: a state which contains a 
superposition of both s- and non-s-wave pairing. For ex- 
ample the s + id and s + d states. In general such a state 
can be written as 

A(r) = A (r) (cos fo>( T )) A s + e w sin foj(r)) A„) (8) 

where r is the quasiparticle lifetime, A(t) is the order 
parameter of the superconductor, Aq(t) gives the mag- 
nitude of the order parameter, A s is a function with a 
magnitude of unity and s-wave symmetry A„ is a func- 
tion with a magnitude of unity and the appropriate non- 
s-wave symmetry and 9 and (p(r) parameterise the su- 
perposition. For clarity we have suppressed all spin and 
momentum labels. We will describe this state as the s + n 
state. 

Naively, it might appear that the s + n state might 
explain the low temperature behaviour of the thermo- 
dynamic and transport properties. If the states had a 
large d-wave component it would appear to have nodes 
at high temperatures, but at low temperatures the small 
fully gapped s-wave part of the order parameter would 
cause an exponential cutoff. However, a more careful 
analysis of the data shows that this scenario is not what 
has been observed, indeed the results of the experiments 
performed to the lowest temperatures suggested nodes in 
the gapi^i 

To describe the effect of disorder on the s + n state 
we will begin my studying the two extreme cases of total 
coherence between the states and zero coherence between 
the states. It will then be seen that all other possibilities 
arc intermediates of these two extremes. 

If there is total coherence between the states, then 
adding disorder does not change the ratio between the 
s-wave and non-s-wave parts of the order parameter, i.e. 
ip is independent of r. It is straightforward to show that, 
subject to this constraint, 
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R 
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where V is the effective pairwise interaction between the 
electrons. For s-wave pairing in the presence of non- 
magnetic impurities^ 



R= 1 



1 
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and one finds that T c — T c0 independent of t, in con- 
firmation of Anderson's theorem. But, for non-s-wave 
pairing R = 1 and we arrive at the AG equation JSJ). 
For an s + n superconductor 



R= 1 



2t\uj„\ 



(11) 



a(ip) is an unknown function, however it is clear that 
a(0) = 1 and a(Tr) = 0. Thus one finds that 
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Thus we find that rigid coherence in an s + n supercon- 
ductor simply 'renormalises' the quasiparticle lifetime in 
the AG equation. 

For a superconductor without coherence between the 
two parts of the order parameter ip varies strongly with 
r and the two parts of the order parameter are in- 
dependent of one another. Thus non-magnetic disor- 
der does not change the bulk critical temperature be- 
cause of the s-wave part of the wavefunction. But non- 
magnetic disorder would reduce the critical temperature 
for the non-s-wave part of the wavefunction. This would 
lead to there being two phase transitions in the pres- 
ence of non-magnetic disorder, the first from the non- 
superconducting state to an s-wave superconductor and 
the second from an s-wave superconductor to an s + n 
superconductor. Two such phase transitions would have 
a clear experimental signature. For example there would 
be two anomalies is the specific heat. This has, to the 
best of our knowledge, never been observed in the layered 
organic superconductors. Therefore we can rule out the 
possibility of s+n superconductivity with zero or, indeed, 
weak coherence between the states on phenomenological 
grounds. 



B. Non-magnetic disorder in other 
superconductors 

The effects of non-magnetic disorder have been care- 
fully observed in several other superconductors. The best 
known case is Sr2Ru04. Mackenzie el alS^ measured T c 
for several samples with varying residual resistivities. As- 
suming the Drude model of conductivity they found the 
variation of T c with po to be in excellent agreement with 
the AG formula. 

Both magnetic (Ni) and intrinsically non-magnetic 
(Zn, Pr, fast electron irradiation) defects lead to the sup- 
pression of T c of YBaCuaOe+a; (YBCO) in line with the 
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AG formulaiSiiSi However, it is known££ that the sub- 
stitution of Zn atoms for Cu atoms in the Cu0 2 planes 
of YBCO can lead to the formation of localised mag- 
netic moments. It is thought that these local moments 
form on the nearest neighbour Cu atoms rather than on 
the Zn site itselfi^l There has been much debateSMS as 
to whether the mechanism for pair breaking in YBCO 
crystals with Zn impurities is local moment scattering or 
potential scattering due to the Zn impurity (of course 
the two mechanisms are not mutually exclusive*^) . Re- 
cent work by Davis et alm&^L indicates that non-magnetic 
scattering is the dominant mechanism by which Zn 
impurities^ lower T c and further that even the magnetic 
impurities (Ni) act primarily as potential scatterers41i 

In the heavy fermion superconductor UPt3 a suppres- 
sion of T c has been observed that is consistent with the 
AG theoryi2ii2£ Surprisingly both magnetic (Ni) impuri- 
ties and non-magnetic (Gd) impurities suppress T c in the 
same way^ In light of the discovery that Ni impurities 
act primarily as potential scatterers in YBCO it seems 
plausible that the same thing may happen in UPt3. Al- 
ternatively some unknown mechanism may be inducing 
local moments around the Gd atoms. This seems un- 
likely as for this to be consistent with the observation 
that Gd and Ni impurities suppress T c in the same way 
this scenario would require the moment induced around 
Gd atoms to be the same as the moment due to Ni atoms. 

The Bechgaard salts, (TMTSF) 2 X where (TMTSF is 
tetramethylteselanafulvalene and X is an anion, for ex- 
ample CIO4 or Re04) are also very sensitive to non- 
magnetic disorder. It has been suggested that this is be- 
cause they are quasi-one-dimensional systems 
Disorder can be induced by x-ray irradiation, alloying 
or by a cooling rate controlled anion disorder transition 
(which we will discuss further below). All of these sources 
of disorder can reduce T c and can even suppress super- 
conductivity altogether and lead to the formation of a 
spin density wavei2&2& 



III. /3-(BEDT-TTF) 2 X 

There are a series of competing ground states in 
both /3-(BEDT-TTF) 2 X and k-(BEDT-TTF) 2 X includ- 
ing antiferromagnetism and superconductivity. By ap- 
plying pressure or changing the anion the ground state 
of these layered organic crystals can be changed, thus it 
is thought that different anions apply different 'chemi- 
cal pressures '^£*A£ For superconducting crystals pressure 
lowers T c . Thus one might expect that by alloying an- 
ions one could observe the same change in T c due to the 
change in 'chemical pressure'. However, if one adds small 
amounts of a second anion the second anion sites will 
act as non-magnetic impurities. Thus, unless the pair- 
ing state is s-wave, alloying anions will suppress T c . The 
suppression of T c should be governed by the AG formula. 

Tokumoto et alm& have produced alloys in the series j3- 
(BEDT-TTF) 2 (I 3 )i- 3; (IBr2) a; . For x = they found that 




R(0)/R(295) 

FIG. 1: The variation of the superconducting transition tem- 
perature of /3-(BEDT-TTF)2lBr2 with the residual resistance 
ratio, R(0)/R(295). The curve is a fit, using the AG for- 
mula assuming the residual resistivity, po cx 1/rt, where rt 
is the quasiparticle lifetime, to the data of Tokumoto et a/i2& 
(squares) who induced disorder by substituting I3 anions for 
IBr2 and Shegolev and YagubskiiSi (circles) who reported re- 
sistivity measurements for several samples. This indicates 
that either both types of impurity induce magnetic moments 
or else the pairing symmetry is non-s-wave. Note that al- 
though we have written R(0) Tokumoto et al. did not actu- 
ally report R(0)/R(295), but R(T C )/R{295) thus their data 
(squares) should be shifted slightly to the left. As Shegolev 
and Yagubskii reported R(T) / R(295) for a range of temper- 
atures near T c we were able to fit to their data to the form 
R(T)/R{295) = R(0)/R(295)+AT 2 and thus determine both 
R(0)/R(295) and T c accurately. 



T c = 7.4 K and for x = 1 they found T c = 2.4 K. Based 
on Anderson's theorem one expects that for s-wave pair- 
ing T c will vary monotonically with x. However, Toku- 
moto et al. found no indications of superconductivity for 
0-2 J; x < 0.7. A natural explanation of this experiment 
is that for small, non-zero values of x the IBr 2 anions act 
as (intrinsically) non-magnetic impurities in /3-(BEDT- 
TTF) 2 l3 and thus quickly reduce T c to zero. Similarly 
for x < 1 the I3 anions act as impurities in /3-(BEDT- 
TTF) 2 IBr 2 and reduce T c to zero for quite small concen- 
trations. This explanation of course requires non-s-wave 
pairing. 

In figure we plot the data for T c against po for j3- 
(BEDT-TTF) 2 (l3)i-x(IBr 2 ) a! with x < 1 from Toku- 
moto et alSSi on the same graph as data for /3-(BEDT- 
TTF)IBr 2 samples^ which have differing residual resis- 
tivities because of impurities accidently induced in the 
fabrication process. The excellent agreement with the 
AG formula is strong evidence against the weak localisa- 
tion scenario. In this fit we assume only that po oc 1/tjv- 
There were not enough data points reported for x < 
to make a similar comparison for /3-(BEDT-TTF)l3. For 
a more detailed discussion of the role of disorder in (3- 
(BEDT-TTF)I 3 see Ref. |H 

It is also interesting to note that the compound 
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/3-(BEDT-TTF) 2 I 2 Br is not superconducting. For 
/?-(BEDT-TTF) 2 X, when X is a trihalide, the three 
positions of the halide atoms are crystallographically 
distinct. In /3-(BEDT-TTF) 2 I 3 the three iodine atoms 
are arranged approximately linearly (which we represent 
by I-I-I) and are clearly indistinguishable particles. In 
/3-(BEDT-TTF) 2 IBr 2 the atoms are arranged Br-I-Br, 
that is to say that the iodine atom is always in one 
particular location. But, in /3-(BEDT-TTF) 2 I 2 Br, 
the atoms can either be arranged I-I-Br or Br-I-I. 
This means that the crystal is intrinsically disordered. 
/3-(BEDT-TTF) 2 I 2 Br is found to have a high residual 
resistivity^ 6 - Thus we propose that it is the intrinsically 
non-magnetic disorder, caused by the two possible 
arrangements of the anion, that suppresses superconduc- 
tivity in /3-(BEDT-TTF) 2 I 2 Br. Further Tokumoto et al. 
observed that no samples with i?(0)/i?(295) > 0.3 from 
any of the alloys /3-(BEDT-TTF) 2 (I 3 )i_ a: (IBr 2 ) K , 
/3-(BEDT-TTF) 2 (IBr 2 ) 1 _ :c (I 2 Br) :c or /3-(BEDT- 
TTF) 2 (I 2 Br)i_ x (l3) a; superconducted. This is exactly 
what one would expect from the AG formalism (c.f. 
figure nj. 

At this stage it may appear that the arguments pre- 
sented above are in contradiction to what is known 
about the cuprate superconductors. These materials 
have d-wave order parameters and yet non-stoichiometric 
compounds often have far higher transition tempera- 
tures than the (stoichiometric) parent compounds (in- 
deed in many cases the parent compound is non- 
superconducting). An excellent example of this is 
La 2 _ a; Sr a ;Cu04 for which optimal doping is x ~ 0.15. 
It was suggested^ that d-wave superconductivity is 
observed in non-stoichiometric compounds because the 
Born approximation is not valid for the cuprates. How- 
ever, it has been showniffi that even in the unitary 
(or resonant) scattering limit which is appropriate for 
the cuprates non-magnetic disorder still destroys d-wave 
pairing in line with the predictions of the AG formula 
and leaves s-wave pairing unaffected. Further unitary 
scattering is the appropriate limiljiSi for the unconven- 
tional superconductor^^ UPt3 and in this material T c 
is suppressed by non-magnetic impurities in a manner 
consistent with the AG formula^ as discussed in section 

HTbI 

However, sofar we have neglected the major difference 
between non-stoichiometric compounds in the organics 
and the cuprates. In the cuprates the change in stoi- 
chiomctry introduces a change in the current carrier con- 
centration. This dramatically alters the ground state 
of the cuprates, this effect is absent in the organics be- 
cause all of the anions have the same electronegativity. 
It should be noted however that, both the cuprates and 
the organics are similarly two dimensional as is attested 
by the ratio of the zero temperature interlayer coher- 
ence length, £j_(0), to the interlayer spacing, a. For 
example in K-(BEDT-TTF) 2 Cu[N(CN) 2 ]Br (Ref. El 
£j_(Q)/a = 5.8/30.016 = 0.19 and in the cuprates^ 
£_l(Q)/& ~ 0.06 — 0.45. Therefore, as both compounds 




defects (%) 

FIG. 2: The variation of the superconducting transition tem- 
perature of /3-(BEDT-TTF)2l3 with the number of impurities. 
The data is taken from Forro et alw^ who induced defects by 
irradiating samples with fast electrons. The curve is a fit to 
the AG formula and equation @. This indicates that either 
the radiation induces magnetic moments or else the pairing 
symmetry is non-s-wave. 

are quasi-two-dimensional and alloying anions suppresses 
T c in /?-(BEDT-TTF) 2 X, it cannot be merely the two di- 
mensional nature of the cuprates which is responsible for 
observation of superconductivity in non-stoichiometric 
compounds. 

It has been shown2& that by alloying anions one can 
introduce enough disorder into the system to suppress 
superconductivity. Assuming that this disorder is non- 
magnetic this rules out s + n superconductivity with any- 
thing other than completely rigid coherence between the 
two states (that is to say that a is independent of ip in 
the language of equation 1)11(1 ). Any other type of coher- 
ence would leave a small residual s-wave component even 
in the presence of very large amounts of disorder. 

Defects can also be induced in materials by irradiat- 
ing them with fast electrons. 86, 10 i Such experiments were 
performed on /3-(BEDT-TTF) 2 I 3 by Forro et al^k who 
noted a marked drop in T c as the number of defects in- 
creased. From figure |21 it can be seen that the fit to the 
AG formula and equation 10} is excellent. Unfortunately 
Forro et al. did not report the residual resistivity of their 
irradiated samples so a comparison with transport the- 
ory cannot be made. Again the excellent fit of the data 
to the AG theory is strong evidence against the weak 
localisation theory. 

We have therefore shown that impurities in /3-(BEDT- 
TTF) 2 X suppress T c via the AG mechanism for three 
sources of impurities: alloying anions, fast electron ir- 
radiation and accidentally created defects from the fab- 
rication process. There is no obvious mechanism for 
any of these methods to form magnetic scattering cen- 
tres. Thus the most natural interpretation is that there 
is non-s-wave pairing in /3-(BEDT-TTF) 2 X and the re- 
duction in T c is due to potential scattering. However, 
there is a strong similarity between the layered organic 
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superconductors and the cupratesi£*i§i£2ii2& in particu- 
lar both are close to an antiferromagnetic phase. As we 
have already noted, the substitution of Zn for Cu in the 
Cu0 2 planes of YBCO leads to the unexplained forma- 
tion of local moments on the Cu atoms neighbouring the 
Zn impurity. Therefore one must consider the possibil- 
ity that a novel mechanism is creating local moments 
in all three of experiments discussed above. This may 
seem unlikely, but until further experimental evidence on 
the nature of the impurities formed in these experiments 
becomes available we cannot use disorder to unambigu- 
ously determine whether or not there is s-wave pairing in 
/3-(BEDT-TTF) 2 X. 

IV. k-(BEDT-TTF) 2 X 

One of the most unusual features of k-(BEDT- 
TTF) 2 Cu[N(CN) 2 ]Br is that T c is dependent on the rate 
at whic h the sample is cooled from T > 80 K (Ref. 

The residual resistivity along the c-axis, po, 
is also dependent on the cooling rate. It would appear 
then, that if one cools K-(BEDT-TTF) 2 Cu[N(CN) 2 ]Br 
quickly one can 'freeze in' disorder, whereas if the cool- 
ing is slower then the disorder can relax out. The ob- 
servation that this disorder suppresses T c implies that if 
the pairing state has s-wave symmetry then the disor- 
der must arise from magnetic impurities, but if another 
pairing symmetry is realised then this disorder may arise 
from non-magnetic impurities. 

There is always a certain amount of intrinsically non- 
magnetic impurities in any given crystal. These 'struc- 
tural' impurities will also contribute to the residual resis- 
tivity, but they only effect T c in the non-s-wave case. We 
denote the quasiparticle lifetime caused by this structural 
disorder by r s . Similarly we will denote the quasiparticle 
lifetime caused by the cooling rate induced disorder by 
t c - 

As non-magnetic impurities do not affect T c for s-wave 
pairing T c is given by JQ) with tm = t c . On the other 
hand both scattering from magnetic and non-magnetic 
impurities contribute to the residual resistivity so we 
might expect 

111 

Po oc — = 1 (13) 

t* t s T c 

where r t is the appropriate quasiparticle lifetime for 
transport experiments. 

The fabrication of different samples will lead to dif- 
ferent values of r s . For s-wave pairing this will cause a 
variation in po but not T c thus one reaches the conclusion 
that different samples cooled at the same rate will have 
different residual resistivities, but the same maximum 
critical temperature. In figure|21we fit the linearised AG 
equation J7J to the data of Su et al*2L We also show the 
effect of varying r s for the s-wave pairing/non-magnetic 
structural impurity scenario, that is that from sample 
to sample the minimum po as a function of cooling rate 
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FIG. 3: Variation of the superconducting transition temper- 
ature, T c , of fi>(BEDT-TTF) 2 Cu[N(CN)2]Br with the inter- 
layer residual resistivity, po. The solid line is a fit to the data 
of Su et alw— (squares). The other lines are predictions of 
the s-wave theory for other samples with different amounts 
of structural disorder and thus a different r s . This structural 
disorder is assumed to be non-magnetic. Thus for s-wave pair- 
ing the structural disorder changes po but does not affect T c . 
The data of Stalcup et alw— then represents a test of the 
theory. It can clearly be seen that the theory does not de- 
scribe the data as both T c and po are changed for all cooling 
rates. This indicates the r s is different for both samples and 
therefore that either the assumption of non-magnetic struc- 
tural disorder is incorrect or the assumption of s-wave pairing 
is incorrect. Note that we have reanalysed the experimental 
data and used a consistent definition of both T c (based on 
when the resistivity falls to half of its normal state value) and 
po (based on a fit to the form p{T) = po + AT 2 ; Matthiessen's 
rule—- was found to be obeyed). 



changes, but the maximum T c does not change. The bro- 
ken lines then show the expected behaviour for different 
samples based on the Su et al. data assuming s-wave pair- 
ing and non-magnetic structural impurities. Also shown 
is equivalent data from experiments performed by Stal- 
cup et al. 108 It is clear that the data from Stalcup et 
al. does not fit with the expectations for s-wave pairing 
and non-magnetic structural impurities. For non-s-wave 
pairing and/or magnetic structural impurities both struc- 
tural disorder and cooling rate induced disorder reduce 
T c . Thus T c is given by © with t n = r t . While the 
residual resistivity is still determined by (|13f) . 

The solid line in figure 0] represents at fit to the data 
of Su et al. The fabrication of different samples will lead 
to different values of r s . This will cause a variation from 
sample to sample in both the minimum value of po and 
the maximum value of T c obtainable by varying the cool- 
ing rate. However as T c and po are both functions of only 
one variable (r t ) the data for all samples will lie on a sin- 
gle line. Thus the broken lines in figure 01 represent the 
prediction of the behaviour of different samples based on 
the Su et al. data assuming non-s-wave pairing and/or 
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FIG. 4: Variation of the superconducting transition temper- 
ature, T c , of K-(BEDT-TTF) 2 Cu[N(CN) 2 ]Br with the inter- 
layer residual resistivity, po- The solid line is a fit to the data 
of Su et aim— (squares). For non-s-wave pairing the struc- 
tural disorder both changes po and T c . The data of Stalcup 
et aim— then represents a test of the theory. The broken lines 
are a prediction of the non-s-wave theory for other samples 
with different levels of structural disorder and thus a different 
t s . It can clearly be seen that the theory describes the data as 
both T c and po are changed for all cooling rates in line with 
the predictions of the AG formula and equation 1221 . The 
dashed portion of the line describes the data of Stalcup et al, 
the dotted line is the prediction for a crystal with even less 
structural disorder. The experimental data and the solid line 
are identical to those shown in figure Note however that 
this figure also represents the prediction for s-wave pairing 
assuming that the structural impurities are solely magnetic 
scatterers. 



magnetic structural impurities. It is clear that the data 
from Stalcup et al. is in excellent agreement with the 
expectations for non-s-wave pairing. 



We stress that this result is based on experiments on 
only two samples. To be conclusive one would require 
the study of many more samples. Further it has been 
arguedAi^ that some measurements of the critical tem- 
perature and residual resistivity in the literatur e) 111 are 
more consistent with the s-wave pairing scenario (figure 
Clearly, a detailed, systematic study is required to 
settle this debate. 

The above work is based on the (reasonable) assump- 
tion that the structural impurities are non-magnetic. 
As we speculated in the case of /?-(BEDT-TTF) 2 X, it 
may be that some novel mechanism of local moment 
formation exists in the layered organic superconduc- 
tors. Applying a hydrostatic pressure or changing the 
anion (X) in k-(BEDT-TTF) 2 X has a dramatic effect 
on the ground state. For example at ambient pressure 
and low temperature k-(BEDT-TTF) 2 Cu[N(CN) 2 ]C1 is 
a Mott-Hubbard antiferromagnetic insulat or. Apply- 
ing a small pressure (~ 200 bar, Ref. EH moves 
k-(BEDT-TTF) 2 Cu[N(CN) 2 ]C1 into an superconducting 
state with properties very similar to those of k-(BEDT- 
TTF) 2 Cu[N(CN) 2 ]Br. Thus it is thought that k-(BEDT- 
TTF) 2 Cu[N(CN) 2 ]Br is close (in anion/pressure space) 
to an antiferromagnetic phase transition^ A possible 
mechanism for the formation of local moments in k- 
(BEDT-TTF) 2 Cu[N(CN) 2 ]Br is that non-magnetic im- 
purities change the local electronic structure by a small 
amount. This small local perturbation could cause the 
formation of a local moment similar to those found in 
the antiferromagnetic phase. A similar suggestion was 
made by Kohno et al—& who considered the competition 
of antiferromagnetic and superconducting ground states 
in Cea;Cu 2 Si 2 with x < 1. In their scenario Ce vacancies 
act as intrinsically non-magnetic impurities, but lead to 
the formation of local moments. At low enough densities 
such magnetic impurities will act as independent, para- 
magnetic spins. As such the impurities' behaviour in a 
magnetic field is governed by the Brillouin function>i2S 



J 



M = N M gvBJz 



1 

2J t 



coth 



1 
2J~ 



2J t I 2k B T 



(14) 



r 



where Nm is the total number of magnetic impurities, 
Ji is the total angular momentum of the impurity and g 
is the usual g-factor. For localised, non-interacting elec- 
trons it is appropriate to take Ji — ^ and g — 2. In which 
case 



M = Nm /-is tanh (j^j ■ 



From (0) we have 
N M 



ZkN{0)\u m \ 2 t c 



(15) 



(16) 



A^(0) is knownii^ because for a quasi- two dimensional 
metal the density of states at the Fermi level is given by 

"(0) = £f (17) 

where m c is the cyclotron mass. In the presence of inter- 
actions Luttinger's theoremii^ for a Fermi liquid ensures 
thaO±i 



N(0) 



(18) 



2vr^ 2 ' 

where m* is the effective mass, regardless of the details 
of the band structure. It is known from Shubnikov-de 
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Haas experiments 1 ^ that, for the (3 or magnetic break- 
down orbit m*/m e — 6.4 and so N(0) = 14.9 eV -1 unit 
cell -1 spin -1 . 

A more difficult problem is estimating um- We can 
make an estimate because of our knowledge of the Mott- 
Hubbard state which is nearby in pressure/ anion space. 
We estimate that um will be of the same order as JV 
where J is the exchange coupling in the Mott-Hubbard 
state and V is the volume occupied by a dimer and an 
anion. This is dimensionally correct and we know that 
in the Mott antiferromagnetic state there is on e spi n per 
dimer. It is estimated that J ~ 40 K (Ref. Ill 7f) and 
hence \um\ = 0.026 eVA 3 . A less theory-laden estimate 
of J can be made from the fact that the Kondo effect is 
not observed in these materials. In the Kondo effect a 
minimum in the resistivity occurs at the Kondo temper- 
ature, Tk, which is given byiA^ 



rr W 1 



(19) 



where W is the bandwidth. For k-(BEDT- 
TTF) 2 Cu[N(CN) 2 ]Br, W = 2{t x + t 2 ) ~ 0.23 eV, where 
ii and i 2 are the nearest neighbour and next nearest 
neighbour hopping integrals respectively? 1 ^ and N(0) is 
given by (|18fl with m* /m e = 6.4. That the Kondo effect 
is not observed implies that Tk < T c < T c q < 12 K from 
the fit in figures [3] and 0] This implies that J < 155 K 



and thus that 



< 0.4 eVA 3 . However, while the 



Kondo temperature is defined for a single impurity, the 
Kondo minimum will not be observable unless there is a 
sufficiently large number of impurities (typically a few 
percent* 1 ^) . 

Substituting 0) into l|16|l we find that 



N M = 



32k f 



3TT 2 hN(0)\u M \ 



;(T c0 -T c ). 



(20) 



For example, Su et alMtl report a maximum variation in 
the critical temperature of T c q — T c = 0.58 K, which leads 
to, as a lower bound (based on J~ 155 K), Nm > 0.03 
impurities per unit cell. For our best guess (J ~ 40 K) 
we find Nm > 0.50 impurities per unit cell. This should 
be sufficient to observe a Kondo minimum and thus the 
Kondo effect places a limit on the number of impurities. 
Substituting J2SJ into (JTBJl we find that, 



32fcf 



M 



fi B 37r 2 /L/V(0)|u A /| 



rO 



T c )tanh(^f). (21) 



Two studies of the variation in magnetisation with 
cooling rate in K-(BEDT-TTF) 2 Cu[N(CN) 2 ]Br have been 
conducted4A2iA2£i Both studies were primarily concerned 
with the weak field limit, but surprisingly even these re- 
sults may tell us something about the presence of mag- 
netic impurities. Taniguchi and Kanoda^S measured 
M(H) at T = 7 K. They found an interesting weak field 
dependence (presumably this is due to vortex dynamics 
as it disappears when the irreversibility line is reached, 



but we will not discuss this here). Above the irreversibil- 
ity line they found that the change in M with cooling rate 
is only weakly dependant on H . (Results were reported 
up to H = 1200 Oe.) Based on the observed cooling 
rate dependence of T c in this sampleiSi we estimate that 
the variation in T c between when the sample is cooled at 
10 K/min and when the sample is cooled at 0.5 K/min is 
0.25 K. This leads to the conclusion that the difference 
in the magnetisation of the two samples due to the mag- 
netic impurities (required in the s-wave scenario) would 
be 1.3 x 10 -4 emu at H=1200 Oe and T=7 K (based 
on our lower bound from the Kondo effect, J = 155 K). 
This is well within the resolution of the experiment (in 
fact this contribution would dominate the observed mag- 
netisation) and is not observed (see figure |5J. Thus the 
experiments of Taniguchi and Kanoda are inconsistent 
with the hypothesis that cooling rate induced disorder 
creates paramagnetic impurities. (However, it is pos- 
sible that paramagnetic impurities are present in the 
sample and that there presence is screened by the su- 
perconducting state.) We therefore suggest that there 
is non-s-wave pairing in K-(BEDT-TTF) 2 Cu[N(CN) 2 ]Br 
and that varying the cooling rate induces non-magnetic 
disorder which causes the variation in both T c and po- 
Again we stress that because the is little data above the 
irreversibility line, H lTl and no data outside the super- 
conducting state, further careful systematic experiments 
are required preferably in the normal state. 

Two groups have investigated anomalies in heat 
capacity 122 ! 123 and thermal expansion 1 ^ at T ~ 8 K 
in k-(BEDT-TTF) 2 X for X = Cu[N(CN) 2 jCl, (Refs. EH 
and [lH Cu[N(CN)?lBr (Refs. CM |l23| and [HI and 
Cu(NCS) 2 (Ref. Il2l . Both groups concluded that 
the anomalies are due to a transition in which disor- 
der becomes frozen into the orientational degrees of free- 
dom in the terminal ethylene groups of the BEDT-TTF 
molecules. This ethylene ordering transition provides a 
natural explanation for the observed cooling rate depen- 
dence of the residual resistivity of k-(BEDT-TTF) 2 X. 
However one should note that such an ethylene order- 
ing transition would result in intrinsically non-magnetic 
impurities and is therefore strong evidence in support of 
our suggestion that the cooling rate induced disorder is 
non-magnetic in nature. 

Terminal ethylene group disorder in k-(BEDT- 
TTF) 2 X is rather similar to the anion disorder observed 
in the Bechgaard salts. In both (TMTSF) 2 C10 4 and 
(TMTSF) 2 Ru04 the anions can occupy two inequivalent 
orientations. Fast cooling leads to partially disordered 
domains, the size of the domains has been shown to be 
proportional to the cooling rate^S^ As mentioned in sec- 
tion [ini] varying the cooling rate can lead to a reduction 
in T c and even the complete suppression of superconduc- 
tivity in favour of a spin density wave. Also note that the 
anion ordering temperature, Tao, is highly dependent on 
which anion is considered. For X = CIO4, Tao ~ 24 K, 
for X = Re0 4 , T AO - 170 K and for X = PF 6 no anion 
ordering transition is observed^ 3 - The nature of the an- 
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FIG. 5: The cooling rate dependence of the magnetisation 
of /t-(BEDT-TTF) 2 Cu[N(CN) 2 ]Br. We plot the difference in 
the magnetisation of the same sample when it is cooled at 
10 K/min and when it is cooled at 0.5 K/min (circles) mea- 
sured by Taniguchi and Kanoda— ■ Also shown is the differ- 
ence in the magnetisation for the same sample when it had 
been annealed at 70 K for 12 hours and when it was cooled 
at 0.5 K/min (diamonds) and the difference in magnetisa- 
tion between when the sample was annealed and when it was 
cooled at 10 K/min (squares). All sets of data were taken at 
T — 7 K < T c . The solid lines are the calculated lower bound 
on the change in the magnetisation at T — 7 K due to param- 
agnetic impurities which produce a 0.25 K change in T c which 
is the estimated change in T c between the sample cooled at 
10 K/min and the sample cooled at 0.5 K/min based on the 
observed cooling rate dependence of this sampled— This lower 
bound is required to ensure the Kondo temperature, Tk < T c 
and thus to be consistent with the fact that the Kondo effect 
is not observed in K>(BEDT-TTF) 2 Cu[N(CN) 2 ]Br. The long 
dashed lines represent the predicted magnetisation assum- 
ing that the interaction energy of the magnetic impurities is 
the same as the observed antiferromagnetic exchange interac- 
tion in the insulating phase of k-(BEDT-TTF) 2 Cu[N(CN) 2 ]C1 
(i.e., J ~ 40 K). The vertical dashed line indicates the ir- 
reversibility line at 7 K, Hi r (T — 7 K), as measured by 
Taniguchi and Kanoda— in the same experiment. Thus we 
see that for H < H- a (T = 7 K) (left of the dashed line) the 
non-trivial vortex dynamics of the system cause a complicated 
variation in the magnetisation, which we do not discuss here. 
However, for H > H ir (T = 7 K) (right of the dashed line) 
the measured difference in the magnetisation is less than that 
required by the Brillouin function. Therefore these measure- 
ments suggest that no paramagnetic impurities are induced by 
varying the cooling rate of this sample. But this conclusion 
requires that the moments are not screened by supercurrents. 



ion order also differs for X = CIO4 and X = ReC>4 (Ref. 
193) . A similar disordering transition is observed— i in the 
organic conductors (DMET) 2 BF 4 and (DMET) 2 C10 4 . 

Of the salts considered here, a variation in T c with 
cooling rate had only been observed in k-(BEDT- 
TTF) 2 Cu[N(CN) 2 ]Br to date. If our hypothesis that the 
variation in T c with cooling rate is due to cooling rate in- 
duced disorder which in turn is due to the ethylene order- 
ing transition in the terminal ethylene groups is correct 
then one would also expect a variation in T c with cooling 



rate in k-(BEDT-TTF) 2 Cu(NCS) 2 as the ethylene order- 
ing transition has been observed in this compoundiiSiiSi 
An ethylene ordering transition has also been observed 
in k-(BEDT-TTF) 2 Cu[N(CN) 2 ]C1 (Ref. How- 
ever, this compound only becomes superconducting un- 
der pressure and it is not known what effect pressure has 
on the disordered ethylene state. Clearly the dependence 
of T c on cooling rate is in need of further investigation. It 
may be of interest to investigate the effect of pressure on 
the ethylene ordering transition, particularly with refer- 
ence to k-(BEDT-TTF) 2 Cu[N(CN) 2 ]C1 and cooling rate 
dependence of the Neel temperature. 

In light of the variation of T c with cooling rate it is 
important that in experiments on the k-(BEDT-TTF) 2 X 
salts the cooling rate is reported regardless of whether or 
not it is varied. Results for T < 80 K lose much of their 
significance if the cooling rate is not known. 

Work by Taniguchi et ali^i^ has raised the pos- 
sibility of inhomogencous phase coexistence between 
antiferromagnetism and superconductivity in dcutcr- 
ated K-(BEDT-TTF) 2 Cu[N(CN) 2 ]Br. There is no ev- 
idence of phase coexistence in fully hydrogenated k- 
(BEDT-TTF) 2 Cu[N(CN) 2 ]Br so phase coexistence can 
be ruled out as the cause of the suppression of T c 
in the hydrogenated compound, which we consider 
here. Further varying the cooling rate of dcutcr- 
ated K-(BEDT-TTF) 2 Cu[N(CN) 2 ]Br and k-(BEDT- 
TTF) 2 Cu[N(CN) 2 ]C1 offers the possibility of studying 
the Mott transition in the presence of disorder with fine 
experimental control over the level of disorder in the sam- 
ple and of varying the level of disorder within a single 
sample. 



V. INTERLAYER TRANSPORT THEORY 

The residual resistivity for interlayer transport in a 
laye red Fermi liquid is given by (see, for example, Ref. 



TTft 4 



1 



Po 



2e 2 m*ct 2 , r t 



(22) 



where c is the interlayer spacing, m* is the effective quasi- 
particle mass and t± is the interlayer hopping integral. 
Thus the assumption that pa oc l/r t l|13|) is justified. 
Substituting into Q we find that 
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(23) 



Thus from our fit to the data of Su et alMH (shown in 
figure |U we have, for K-(BEDT-TTF) 2 Cu[N(CN) 2 ]Br, 



T,, 



and 



11.7 K 



0.9 Ocm. 



(24) 



(25) 
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Taking m* = 6.4m e (Refs. EM and ITT1 and c = 
30.016 A (Ref. H3) we have t± = 0.022 meV. How- 
ever, we note that m* was determined for the 13- 
sheet (that is the magnetic breakdown orbit) only 
whereas here we are considering an effective one band 
model. Nevertheless this value is in excellent agree- 
ment with an independent determination of t± from 
angular-dependent magnetoresistance (AMRO) experi- 
ments. Although t± has not been measured experimen- 
tally in K-(BEDT-TTF) 2 Cu[N(CN^]Br, for k-(BEDT- 
TTF) 2 I 3 t ± « 0.016 meV (Ref. [Hi and for k-(BEDT- 
TTF) 2 Cu(NCS) 2 t ± w 0.04 meV (Ref. Hal . 

For /3-(BEDT-TTF) 2 IBr 2 (see figure HJ we find that 
T c o = 3.0 K. Tokumoto et a/. 96 reported that the room 
temperature resistivity of their samples was p(295) = 
(5.0 ± 2.5) x 10~ 2 Slcm. Therefore 



e 2 m*ct]_ 
Ak B h 3 



40 ± 20 Slcm. 



(26) 



Taking m* = 4.2m e (Ref. Q and c = 15.291 A (Ref. 
0) we have t± = 0.26 ± 0.07 meV. Note that this 
is an order of magnitude larger than for k-(BEDT- 
TTF) 2 Cu[N(CN) 2 ]Br. However, this value is also in 
agreement with previous estimates from de Haas-van 
Alphcn experiments. Wosnitza et alM^ showed that for 
/?-(BEDT-TTF) 2 IBr 2 , t±/E F « 1/280 they also found 
that k F ~ 3.46 x 10 9 m _1 . Therefore taking m* = 4.2m e 
(Ref. Il4|) again and assuming a cylindrical Fermi surface 



h 2 k 2 F 

2m* 



(27) 



one finds that t± ~ 0.35 meV in excellent agreement with 
our result. 

The agreement between t± calculated from our fits 
via equation (|22fl and the values found from AMRO ex- 
periments for both /3-(BEDT-TTF) 2 IBr 2 and k-(BEDT- 
TTF) 2 Cu[N(CN) 2 ]Br is further evidence that in these 
compounds T c is suppressed by the AG mechanism and 
not by weak localisation. 

It has recently been shownS 8 . that the ob- 
served variation of T c and po for alloy /3-(BEDT- 
TTF) 2 (l3)i-x(IBr 2 ) a; for small x predicted by JHI is 
consistent with the observations of Tokumoto et a/. 96 
Note that this theory has no free parameters once the 
T c0 (this work and Forro et al. 105 ) and t± (AMRO 
experiments^^,) have been determined. 

The agreement between our calculated values of t± 
and those measured in AMRO experiments indicates that 
if there is an s + n state then the s-wave component 
(cos(y))A s ) is small (see section iHl . (Or more strictly 
that a is small, c.f. equation (|12|) .) It therefore appears 
unlikely that the layered organics are s + n superconduc- 
tors. 



VI. DISCUSSION 

This study of the effects of disorder on the layered or- 
ganic crystals /3-(BEDT-TTF) 2 X and k-(BEDT-TTF) 2 X 
has shown that disorder has the potential to differenti- 
ate between s-wave and non-s-wave pairing states. But, 
more experiments are needed. This is largely because 
none of the experiments that we have discussed in this 
paper were designed to study the pairing symmetry. In 
this section we will explore what the unresolved issues 
are and how they could be resolved. 



A. Sample variation 

Perhaps the simplest test for unconventional super- 
conductivity is to study the variations in the supercon- 
ducting critical temperature reported in the literature. 
Crystal growers go to great lengths to avoid the inclu- 
sion of magnetic impurities, but the inclusion of non- 
magnetic impurities 178 is harder to avoid. For exam- 
ple the first reports of superconductivity in Sr 2 Ru04, 
which is widely considered to have an unconventional 
(triplet) pairing symmetry, indicated that T c — 0.93 K 
(Ref. [135J). However, sample quality was rapidly im- 
proved and it is now believed that the maximum critical 
temperature, T c o = 1.5 K (Ref. I&ij) has been achieved. 
Thus, for Sr 2 Ru04, T c has increased by over 50% since 
the first report of superconductivity. In contrast, con- 
sider MgB 2 . The first reporliiSi of superconductivity 
quoted T c = 39 K. No significant increase in T c has 
been reported thus far. This is evidence for s-wave pair- 
ing in MgB 2 . Further, doping MgB 2 with U does not 
significantly alter T n ( doping with 1 wt% U reduces T c 
by < 0.5%, Ref. Il37|) . This is in agreement with the 
emerging consensus that MgB 2 is a strong coupling s- 
wave superconductoriiS 8 . (For a fuller discussion of the 
effects of disorder in MgB 2 see Ref. [l39j.) 

The initial reports of superconductivity in k-(BEDT- 
TTF) 2 Cu[N(CN) 2 ]Br quote T c = 10.8 K (Ref. HM . 
While we have shown that T cQ = 11.7 K. k-(BEDT- 
TTF) 2 Cu(NCS) 2 also shows wide variation in T c from 
sample to sam ple s ome authors have reported T c as low 
as 8.7 K (Ref. Il4l|), wh ile other studies have found that 
T c = 9.3 K (Ref. fill . One complication arises from 
the variety of definitions used to determine T c . Tak- 
ing a resistivity measurement as an example, the T c can 
be defined in a variety of ways: (i) the temperature at 
which p first begins to deviate from the Fermi liquid 
form (p(T) = po + AT 2 ), (ii) the highest temperature 
at which p(T) = 0, or (iii) the midpoint of the tran- 
sition, i.e., the temperature at which the p(T) is 50% 
of the Fermi liquid value. For example, definitions (i) 
and (ii) give a difference of ~ 1 K for the data reported 
by Stalcup et al£& about the value T c = 11.6 K (de- 
fined by method (iii), which we use throughout this pa- 
per). The large variations in T c noted above (8% for 
K-(BEDT-TTF) 2 Cu[N(CN) 2 ]Br and 7% for k-(BEDT- 
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TTF) 2 Cu(NCS) 2 ) are probably too large to be explained 
by subtle variations in the definition of T c and are there- 
fore unlikely to occur for s-wave pairing although this is 
far from conclusive. 

/?-(BEDT-TTF) 2 I 3 shows a strong variation in T c . In 
the j3n phased Kahlich et alm^ reported that T c varied 
between 4.5 K and 7 K depending on which sample they 
measured. This represents a 36% variation in T c . This is 
also suggestive of non-s-wave pairing. 

The wide variation in T c from sample to sample is 
something that great care should be taken over in ex- 
periments designed to study the isotope effect. In par- 
ticular any such experiments need to demonstrate that 
crystals that are nominally identical do indeed have a 
highly reproducible T c . If this is not possible then the 
T c variation within nominally identical samples needs to 
be carefully accounted for. For example, by studying 
the sample dependence of the residual resistivity across 
a range of nominally identical samples and using this to 
calibrate the impurity dependence of the various isotopes. 



B. Measurement of the scattering time 

Disorder would be a much more powerful probe if there 
existed a method by which the scattering time could be 
measured directly. The most obvious techniques for this 
are Shubnikov-de Haas and de Haas-van Alphen exper- 
iments. These quantum oscillation experiments measure 
the quasiparticle lifetime via the Dingle temperature, Tp. 
However, the lifetime determined by quantum oscillation 
experime nts, t„, i s not the same as the transport lifetime, 
Tt (Refs. Il43l and ll44tl . Even in the best experiments on 
elemental metals, it is not at all unusu al for Tt to be 10 
or even 100 times larger than r q (Ref. Il45|) . In partic- 
ular Td and hence T q are known to be very sensitive to 
the mechanical state of the sample. A slight deformation 
caused by, for example, handling the sample can lead to 
dramatic increase in Td (decrease in r q ), whilst hardly 
affecting the electrical resistivity (pg oc 1/tj). Given the 
large compressibility of the layered organic superconduc- 
tors T q is unlikely to be the same as r t . 

In its immediate location a dislocation acts just like a 
line of point defects and thus contributes equally to both 
transport and quantum oscillation experiments. How- 
ever, the long-range strain field produced by a dislocation 
only produces very small angle scattering (as the electron 
wavelength is smaller than the characteristic length scale 
of the dislocations) . Therefore the long-range strain field 
contributes negligibly to the transport lifetime but can 
strongly suppress r q even at relatively low dislocation 
densities. 

A sample which is nominally a single crystal is in fact 
made up of a large number of grains. One can think of 
this mosaic structure of grains as a certain pattern of 
dislocations. In this way it is clear that mosaic structure 
causes highly anisotropic scattering and thus leads to the 
suppression of r q . 



TABLE I: Comparison of the transport/ Abrikosov-Gorkov, 
cyclotron resonance and quantum oscillation quasiparticle 
lifetimes (t ag , r cr and r q respectively). As r is clearly a 
highly sample dependent property this table is not intended 
to report universal results but is indicative of general trends. 
NS indicates a non-superconducting compound for which t ag 
cannot be determined. The reported t ag is based on the sam- 
ples used for the experiments discussed in this paper (or in 
Ref. [l4]] in the case of ZrZn 2 ). We have abbreviated BEDT- 
TTF to ET iu this table. 



Material 


r t AG (ps) 


Tct (ps) 


Tq (ps) 


K-(ET) 2 Cu[N(CN) 2 ]Br 


2.5-20 


? 


0.5-0.6i2& 


/3-(ET) 2 IBr 2 


0.6-1.8 


? 




/3-(ET) 2 I 3 


? 


? 


24 M8 


0-(ET) a Ia 


? 


15-36 149 


0.6-1.5i2a 


a-(ET) 2 KHg(NCS) 4 


?° 


15 152 


0.3-0.5±S 


a-(ET) 2 NH 4 Hg(NCS) 4 


? 




^IM 


(ET) 2 Br(DIA) 


NS 


4.6-5.5^ 


1 


(ET) 3 C1(DFBIB) 


NS 


5.6^ 


1 


Sr 2 Ru0 4 


6-38' 


10-40±S 




ZrZn 2 


~6ltt 


? 


0.3±S 



a a-(ET) 2 KHg(NCS) 4 is only superconducting under pressuraASi 
''The sample measured by Hill et aU& had T c = 1.44 K (for which 
the AG formula gives r t AG = 37.9 ps based on T c q = 1.52 K, the 
value found from fitting the AG formula to the data of Mackenzie 
et a/™). r t AG = 6.25 ps based on T c = 1 K, the value reported in 
the de Haas-van Alphen experiments™™ 

Many previous authors have pointed out the difference 
in the transport and quantum lifetimes. However, Hill 146 
noted a similar difference between the lifetime observed 
in cyclotron resonance experiments, r cr , and the quan- 
tum lifetime. It is therefore interesting to compare the 
lifetimes from cyclotron resonance and quantum oscilla- 
tion experiments with the transport lifetime determined 
from the linearised AG equation Q and the value of T c q 
found from the fit to experiment, r t AG , (see Table [|J. 

We see that T t AG ~ T cr across a broad range of (BEDT- 
TTF) 2 X salts, while T q is consistently an order of mag- 
nitude smaller. This suggests that scattering events are 
not the dominant contribution to Td (c.f., Singleton et 
ali^). It presents the intriguing possibility that cy- 
clotron resonance experiments could be used to probe 
the quasiparticle lifetime and thus directly compare the 
experimental T c with the predictions of the AG equa- 
tion. Indeed cyclotron resonances have already been 
observedi££*i££ in Sr 2 Ru0 4 . The observed cyclotron res- 
onance lifetime is larger than the observed lifetime in de 
Haas-van Alphen experiments, but this may be partly 
explained by the much higher T c of the sample used for 
the cyclotron resonance experiments. Excellent agree- 
ment is found between the measured cyclotron lifetime 
and the lifetime calculated from the AG formula. Clearly, 
a systematic study of how the cyclotron resonance life- 
time (and indeed the quantum oscillation lifetime) varies 
with T c is needed. Sr 2 Ru0 4 would be an ideal mate- 
rial for such experiments as the AG formula is seen to 
be obeyed^ and good quality quantum oscillation^ 
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and cyclotron resonance experiments^^ can be per- 
formed. Alternatively the AG behaviour of k-(BEDT- 
TTF) 2 Cu[N(CN) 2 ]Br would make it an excellent mate- 
rial for such a experiment. This is particularly elegant 
as the cooling rate can be used to vary the disorder and 
hence the scattering lifetime, so the experiment could be 
performed on a single sample. Measurements of the vari- 
ation of the Dingle temperature with cooling rate have 
already been made>i2& 

Kartsovnik, Grigoriev and coworkersi^LiSi have also 
investigated the relationship between the quasiparticle 
lifetimes caused by solely microscopic scattering events, 
and the lifetime extracted from the Dingle temperature 
which also contains the effects of macroscopic inhomo- 
geneities. They have shown that the slow oscillations 
observed in quantum oscillation experiments on quasi- 
two dimensional metals are damped by a modified Dingle 
temperature, Tjj, which is not effected by macroscopic in- 
homogeneities. For experiments performed on /3-(BEDT- 
TTF) 2 IBr 2 they found an order of magnitude difference 
between r q (1.5 ps) and the lifetime derived from T£,, t* 
(8.1 ps). 

The Fermi velocity, Vp, for both the (3 and k poly- 
morphs is typically vp ~ 10 5 ms _1 (see section Ivl and 
Ref. Il6l[) . And we have shown here that a quasiparticle 
lifetime of the order Tj ~ 0.1 ps is required to com- 
pletely suppress superconductivity. Thus the mean free 
path, I = VpTt, is typically I > 10 nm (c.f. Ref. 126^ . 
The interlayer coherence length, £y , is typically a few nm 
(c.f. Ref. yjj). Thus these materials are in the clean limit 
even when superconductivity is completely suppressed by 
disorder. This is further confirmation that the AG mech- 
anism is responsible for the suppression of superconduc- 
tivity in these materials. 

We will conclude this section by outlining a series of 
experiments that could determine if the disorder in the 
layered organic superconductors is due to scattering from 
localised moments or potential scattering. These exper- 
iments therefore have the potential to rule out s-wave 
pairing. 



C. Identification of the pairing symmetry 

Comparatively little attention has been focused on the 
pairing symmetry of /3-(BEDT-TTF) 2 X so we will start 
by considering this crystal structure. All of the meth- 
ods of creating disorder considered in this paper (namely 
fast electron irradiation, alloying anions and accidental 
disorder from the fabrication process) should be revis- 
ited and studied in more depth. Both figures ^ an( l [21 
need more data points. Therefore it is most important 
that the entire AG is mapped out. In particular it is 
important to observe the complete suppression of super- 
conductivity by very small amounts of disorder that is 
a unique feature of the AG formalism. Careful observa- 
tion of the entire AG curve is required to rule out other 
mechanisms for the suppression of T c such as weak lo- 



calisation, interband scattering, changes in the pairing 
interaction or the macroscopic coexistence of supercon- 
ducting and non-superconducting phases. All of these 
mechanism for the suppression of T c produce markedly 
different relationships between T c and po and thus would 
be ruled out by the observation of the entire AG curve 
and in particular the complete suppression of T c by mod- 
erate amounts of disorder which is not caused by any of 
the other mechanisms for T c suppression. Forro et alMl 
did not measure the resistivity of their irradiated sam- 
ples. It is important to know the resistivity for several 
reasons: (i) it allows for easy comparison with other tech- 
niques, in particular it allows a consistent definition of T c 
to be used, (ii) it provides a check on the estimation of 
the number of defects produced and (iii) it allows for the 
calculation of t± and thus for a further check that AG 
theory is indeed relevant. All of these methods should 
also be applied to k-(BEDT-TTF) 2 X. 

The next step is to discover whether any of the meth- 
ods for producing impurities create magnetic scatterers. 
One way to do this is to measure the magnetisation 
as Taniguchi and KanodaiSS have for cooling rate in- 
duced disorder in K-(BEDT-TTF) 2 Cu[N(CN) 2 ]Br. How- 
ever, this experiment should be repeated in the metallic 
state. This suggests that paramagnetic impurities are 
not induced by varying the cooling rate of k-(BEDT- 
TTF) 2 Cu[N(CN) 2 ]Br. 

Here we will consider alternative experiments which 
could be used to search for magnetic impurities. We 
will describe these experiments in the context of cooling 
rate induced disorder in K-(BEDT-TTF) 2 Cu[N(CN) 2 ]Br. 
However, the generalisation of these experiments to the 
other methods of producing disorder is straightforward. 
Cooling rate induced disorder experiments are particu- 
larly elegant as the level of disorder can be controlled 
within a single sample. This reduces systematic errors, 
for example, by far the largest source of error in mea- 
suring po comes from measuring the samples dimensions, 
such errors cancel in cooling rate induced disorder exper- 
iments. 

Muon spin relaxation (^tSR) experiments are capable 
of detecting localised magnetic moments. 162 If local mo- 
ments are produced, then the muon spin relaxation rate 
would vary as a function of cooling rate. Clearly the 
muon spin relaxation rate is changed by the supercon- 
ducting state. As T c and presumably H C 2 are changed by 
the cooling rate it is important that these experiments be 
done in the non-superconducting state, either above T c 
or above if C 2- As the ethylene ordering transition occurs 
at T ~ 80 K and T c ~ 10 K any local moments should 
be well formed several Kelvin above T c . 

Nuclear quadrupole resonance (NQR) experiments 
have been used to observe the formation of local mo- 
ments in La 2 _ a; Sr a; Cu04 (LSCO) for x = 0.06 (Ref. 
Il63j) . As perviously discussed NMR measurements 
have observed localised moments induced by Zn im- 
purities in YBCO.- 8 - Therefore studying the change in 
1/Ti with cooling rate in K-(BEDT-TTF) 2 Cu[N(CN) 2 ]Br 
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could determine whether or not local moments are 
formed. The change in 1/T\ as a function of cooling 
rate has been measured in 98% deuterated k-(BEDT- 
TTF) 2 Cu[N(CN) 2 ]Br. No change in l/T x was observed 
until below 30 K, in particular \jT\ is independent 
of cooling rate near 80 K where the ethylene ordering 
transition occursii 6 ^ However, fast cooling of deuterated 
rc-(BEDT-TTF) 2 Cu[N(CN) 2 ]Br drives the ground state 
from superconductivity to an antiferromagnetic Mott 
insulatoi»i22i±2& (which causes the observed difference in 
1/Ti below 30 K). Therefore this observation does not 
rule out the possibility of local moments in hydrogenated 
K-(BEDT-TTF) 2 Cu[N(CN) 2 ]Br. Wang et al. 165 carried 
out an electron spin resonance (ESR) experiment on k- 
(BEDT-TTF) 2 Cu(NCS) 2 . Wang et al. saw no signal at- 
tributable to Cu(II) species at any temperature although 
they do not comment on other sources of magnetic im- 
purities. Therefore it is reasonable to hope that further 
ESR studies may shed some light on the issue of magnetic 
impurities. 

The techniques, outlined here, for using intrinsically 
non-magnetic disorder to probe the superconducting 
state are clearly more general than the context of [3- 
(BEDT-TTF) 2 X and k-(BEDT-TTF) 2 X that we have 
examined here. Disorder has already been used to study 
Sr 2 Ru04 (Ref. I&ih (although we should note that no 
experiments have been performed to rule out magnetic 
impurity formation in this material). Similar results 
for UPt3 (Ref. l92|) appear to have gone largely un- 
noticed. Clearly more careful analysis of this work is 
required. These methods could also be extended to 
other heavy fermion superconductors. There are sev- 
eral other quasi-two-dimensional organic superconduc- 
tors (such as A-(BETS) 2 X, 6»-(BEDT-TTF) 2 X and 0"- 
(BEDT-TTF) 2 X) which may be unconventional super- 
conductors. Disorder would appear to be a powerful tool 
for the investigation of the superconducting state in these 
materials. 

But, the study of disorder, perhaps, is most power- 
ful when used to identify s-wave pairing. An excellent 
example from the recent past is the high temperature 
superconductor, MgB 2 , which appears to be a phonon 
mediated s-wave superconductor This could be con- 
firmed by careful study of the effects of disorder and 
showing that disorder can b e int roduced with only a 
small change in T c (c.f. Ref. Il39f) . This could also be 
applied to other superconductors which are suspected of 
being s-wave, in particular superconductors suspected of 
having anisotropic s-wave order parameters, such as the 
borocarbidesiiSi 



VII. CONCLUSIONS 

We have considered the effect of impurities and dis- 
order on the superconducting critical temperature in 
/3-(BEDT-TTF) 2 X and k-(BEDT-TTF) 2 X. We have 
shown that various sources of disorder (alloying anions, 96 



fast electron irradiationjiS^ disorder accidentally pro- 
duced during fabrication^ and cooling rate induced 
disorde»iS£*i2£) lead to a suppression of T c that is well 
described by the Abrikosov-Gorkov formula. This is con- 
firmed not only by the excellent fit to a theory with only 
two free parameters, but also by the excellent agreement 
between the value of the interlayer hopping integral, t±, 
calculated from this fit and the value of t± found from 
AMRO experiments. This makes a pairing state with a 
superposition of s-wave and non-s-wave components ex- 
tremely unlikely Although such and s + n state cannot 
be strictly ruled out, the s-wave part of the wavefunc- 
tion must be very small and the coherence between the 
s-wave and non-s-wave parts of the wavefunction must 
be completely rigid (a(ip(r)) = a€l). The agreement 
between the measured and calculated values of t± effec- 
tively leaves T c q as the only free parameter in the theory. 
In practice one has very little choice over the value of 
T c o, so the agreement with experiment is found from an 
essentially parameter free theory. The AG formula de- 
scribes the suppression of T c by magnetic impurities in 
singlet superconductors, including s-wave superconduc- 
tors. However T c is suppressed in exactly the same way 
by non-magnetic impurities in a non-s-wave supercon- 
ductor. We therefore have shown that there are only 
two scenarios consistent with the current state of experi- 
mental knowledge. We summarise these scenarios below. 
The task is now to discover whether the impurities are 
magnetic or non-magnetic. 

Scenario 1: d-wave pairing. If the disorder induced 
by all of the four methods considered in this paper is, as 
seems most likely, non-magnetic then the pairing state 
cannot be s-wave. Triplet pairing is ruled out by the 
combination of the three experiments discussed in sec- 
tion [i] 20 ' 54 > 57 Therefore we know that the angular mo- 
mentum, I, of the Cooper pairs is even. If the disorder 
does turn out to be non-magnetic then this implies that 
I > 2. In which case Occam's razor suggests that d- 
wave pairing is realised in both /3-(BEDT-TTF) 2 X and 
k-(BEDT-TTF) 2 X. 

Scenario 2: a novel mechanism for the forma- 
tion of local magnetic moments. Given the proxim- 
ity of /3-(BEDT-TTF) 2 X and k-(BEDT-TTF) 2 X it the 
Mott-Hubbard antiferromagnetic state in anion/pressure 
space, it is possible that disorder induces local magnetic 
moments. Further Taniguchi et a/piSLiS have suggested 
that varying the cooling rate can to the macroscopic co- 
existence of superconductivity in deuterated k-(BEDT- 
TTF) 2 Cu[N(CN) 2 ]Br. Although there is no evidence for 
anything but a spatially uniform superconducting state 
in the hydrogenated compoundfi2Li2S which we have 
considered here, these experiments would not detect iso- 
lated magnetic impurities. On the other hand we have 
shown here that the work of Taniguchi and Kanoda— - 
is inconsistent with the theory that disorder modulates 
the local electronic structure and thus moves single sites 
or small clusters of sites into a state, analogous to the 
Mott-Hubbard insulating state with localised electrons, 
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which can act as magnetic point scatterers. However, 
only a little data was reported in the relevant magnetic 
field range so further work is needed to rigourously test 
this scenario. 

We have suggested experiments to differentiate be- 
tween these scenarios. Such experiments would either 
discover a novel mechanism for the production of lo- 
calised magnetic moments or determine that the su- 
perconducting order parameter is d-wave in /3-(BEDT- 
TTF) 2 X and k-(BEDT-TTF) 2 X. 
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what we term a non-s-wave is that the phase of the order 
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this can be the case for extended s-wave states on an 
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in T c is predicted, and indeed observed, for s-wave su- 
perconductors. The two main points to note about this 
change in T c are: (i) the number of impurities required 
to produce a small change in T c in an s-wave supercon- 
ductor would be enough to completely suppress super- 
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